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ON THE ASYMPTOTIC CHARACTER OF THE SOLUTIONS OF 

CERTAIN LINEAR DIFFERENTIAL EQUATIONS 

CONTAINING A PARAMETER* 

BY 

GEORGE D. BIRKHOFF 

The aim of the present paper is to develop the asymptotic character of the 
solutions of linear differential equations of the form 

for large \p\. The functions a { (x, p) are analytic in the complex parameter p 
at p = oo and have derivatives of all orders in the real variable x. Schlesinger f 
has proved asymptotic properties for lim p = oo on some fixed ray arg p = a; 
in this paper we prove similar properties for a region = arg p = yjr, but by a 
different method. In 1837 Liouville J treated the important special case 

^+ \_p 2 + g(x)]* = o 

when p is real, the first problem of the kind to be considered ; the method of 
attack used in this paper is of a similar nature. 

It is purposed to make an application of the results here obtained to boundary 
value and expansion problems in a second paper. 

I desire to make acknowledgment of the kind encouragement and valuable 
suggestions received from Professor E. H. Moore, for both papers. 

We consider functions z(x, p) of a real variable x on the closed interval 



*Part of a paper presented to the Society (Chicago), Mar. 30, 1907, under a different title. 
Received for publication October 23, 1907. 

t Mathematieche Annalen, vol. 63 (1907), pp. 277-300. This paper appeared after the 
writing of the present paper. See also a paper by J. Horn, Mathematische Annalen, vol. 
52 (1899), pp. 271-292. 

JLiouville's Journal, vol. 2, p. 16. Cf. in particular \ 3. 
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(a,b) and of a complex parameter p, \p\ > R, which satisfy a linear differ- 
ential equation of the form 

W di +pa ^ X,P ">d^ x+ ■•■+/ , "«o( a! >P) 2i = - 

Here we assume that the functions a t (x, p) are bounded by the inequalities 

(2) \a.(x,p)\SM (a^x^b,\ P \>B), 
and that 

(3) a.(x,p) = ±a. J (x)p-'. (\p\>B). 

The coefficients a., of this last series we assume to be continuous with continuous 

v 

derivatives of all orders on ( a , b ) , but not necessarily real. As a consequence 
of (2) and (3) the a.(x, p) are continuous in x andp, for the series (3) are series 
of continuous functions of x and p converging uniformly for \p\ = H > M. 
Further we postulate that the n roots 

w x (x), w 2 (x), •••, w n (x) 
of the equation 

(4) w" + a n _ 10 (x)w«- 1 +■■■ + a w (x) = 

are distinct for every x on (a, b). 

By fi(p>) we denote the real part of the complex number fi. 

Definition. By a region S of the p-plane we understand a region for 
which the indices 1 to n can be so arranged that 

(5) #!>,(»)] ^ -»!>,(»)] ^ • • • = * O.OO] 

for every x on (a,b) and every p on S. 

If p = p Q is an S point, so that (5) obtains for p = p a , it is clear that (5) 
holds if only arg p = arg p . By virtue of this fact the half line arg p = arg p 
belongs to the same region S. For a given x the relations (5) define a certain 
closed sector, 

x ^axgp^f x , 

containing the half line arg/j = arg/o . The largest closed sector common to 

all these sectors, 

(6) flSiargp^, 

forms the region S, degenerating to a half line when = i}r. The funda- 
mental theorem refers to the solutions of (1) on these regions S, when such 
regions exist. 
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An important set of inequalities on a region S is 

(7) I /J? wmiU < I e <\/N(0* <...<] e ofi w ^ dt I ( a < a < ^ < 6 ) . 

In an expression of the form 

«"«*> [«„(«>"* + ^(a;)^- 1 + • • •] 
the term 

e {x)p- h 

shall be called the principal term. 
Finally we introduce the notation 

4>^(x,p) = P - k ^~ k 4>(x, P ), 

so that (1) becomes 

(8) A(z) = z« + a n _ Y {x, p)z^ + ■■■ + a (x, P )z=0; 

and also the notation E as a generic notation for functions of p and other vari- 
ables bounded for \p\ = i?° , when _fi?° is sufficiently large. 

Lemma I. For every value of ifrom 1 to n inclusive there exist an infinite 
number of functions u m (x) , u a (x), • - • , continuous and with continuous deriv- 
atives of all orders, such that u i0 (x) does not vanish at any point of (a, b) 
and such that if the functions 

u i {x,p) = e Ja X,u(x)p-J 

3=0 

be substituted for z in the expression A(z), the coefficient of 

e p v.=o,i, ■■•,«; 

in the expression thus obtained vanishes identically. 
To prove this, write 

/ \ pf*MOdt 

u(x,p) = e- /a v(x,p), 
where w(x) is some definite one of the n roots w { (x). We find then 

uV\x,p)=e 9 f* wmt \_w{x)v(x, p) +« [1] («, p)} , 

«P](a>, p) = eS"~ lt *[([u>(x)y+-^v>(x)\v(x, p)+ 2w(x)v^(x, p)+v™(x, p)\ , 
and, in general, 

p /~ x wit)dt3 , „, 

(9) uf^ar.pjse^" £«•*(*, /»)t>™(*, p), 
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in which a. 4 (a5, p) is a polynomial in 1/p of degree j — 1 at most, whose 
coefficients are functions of x continuous on (a, 6) and with continuous deriva- 
tives of all orders. We note that the principal term in the coefficient of v(x, p) 
in (9) is 

(10) «„ ( x )= ["(»)]'. 

and the principal term in the coefficient of v w (x, p) is 

(11) ^WIX*)?'- 1 ; 
also we note that 

a..(x,p) = l. 
Hence we conclude 

(12) A[u(x,p)]=e'' j7w(t) ' U I [v(x,p)] 
where 

(13) l(z) = ssM + 5._ 1 (aj, /»)«[— »1 + • • • + a (x, p)z, 

the a. (as, /j) being conditioned as the a.(x, p) are in equations (2) and (3). 
For convenience we write 

«n( a5 '/ , )= 1 > « ; ( a5 »P) =0 (*>«)• 

We see from (10), (11) that if we place 

a t (x,p) = Y,%.(x)p~ J , 
then 

(14) a 09 (x) = [w(x)Y + a n _ li0 (x) [w(x)y~ l +■■■+ a m (x) = 

and 

(15) a 10 (x) = n [w(x)y~ l + (w - 1)^ (x)[>(x)]"- 2 + • • • + a l0 (x), 

so that a 10 (x) =j= on (a, b), the w roots ^(x) being distinct. 
If then we write in (12) 



v(x, p)= ^u.(x)p J , 

J=0 



we find the condition that the coefficient of 



in A [«(cc, p)] vanishes to be 

(16) £ a (x)- f - l u k {x) = 0. 



e Ja P~° (s = 0, 1, •••,«) 



J+fc+fcS 



The equation (16) is true for s = by (14). For s=)=0 we can write (16) in 
the form 

£ k<s-l gl 

(17) a u (x)-^u s _ l {x) + a (ix {x)u t _ l {x) + T, = %.(x)^ xl u k (x) = 0, 
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since, if in (16) k = s, then I =j =. and the term corresponding to this set 
of values has a coefficient a oo (x) = 0; if k = s — 1, we have either 1 = 1, 
j = or I = , j = 1 , and the corresponding terms are the first terms of (17). 
It appears then that w a _,(a;) can be determined in terms of 
u i-t(. x )> u ,-z( x )i "'■> u o( x ) as a solution of a certain linear differential equa- 
tion of the first order which has no singular points on (a, b) since a l0 (x) 4= 0. 
Thus u (x), u Y (x), •••, u m _ l {x) are obtained in succession from (16) for 
s = 1 , 2, • • •, m. For each w { (x) we obtain in this way a sequence of func- 
tions u i0 (x), «.,(&), • • ., such that if the expressions 

u i (x,p) = e Ja 2^ u i,( x )P 

be substituted for a in -4(a), the coefficient of 

e /a p~' (i = l, 2, •••, n; s = 0, 1, ••-, m) 

vanishes by (17) since the conditions (16) are now satisfied for s = l,2,---,m. 
Furthermore the differential equation for w i0 (se) is homogeneous, so that by 
taking for u m (x) a solution which is different from zero at one point of 
(a, b), we are sure that u i0 (x) does not vanish at any point of (a, b). 
Since the «•,(*) were continuous with all of their derivatives, the functions 
m. (x), ti a (x), ■ ■ • are also continuous with all of their derivatives. The sequence 
of functions u i0 (x), u n (x), • • • has then the properties stated in the lemma. 
If the formal developments 

P /I Wi(t)dt ^— * . . 

« E«^)r' (i=i, 2, •••,«) 

converged and admitted of ?i-fold term by term differentiation, we should have 
in them the asymptotic developments of solutions of our differential equation 
which we desire. This however will not in general be the case. We shall, 
nevertheless, show in the theorem below that by breaking off after m terms of 
these series (i. e., by retaining the part u^x, p)) and putting in certain remain- 
der terms we can get true solutions of our differential equation. In order to 
prove this we must first prove another lemma in which the form of the dif- 
ferential equation is established which is satisfied by the sums of these m terms 
(i. e., u^x, p)) without any remainder : 

Lemma II. The homogeneous linear differential equation of order n with 
n solutions « t .(x, p) has the form 

B(z) = Z W + &,_,(£, p)z^ + ■■■ + b (x, P )z = 0, 

where for | p | > i?° 

(18) \b { (x,p)\^M°; h i {x,p) = f^b i .(x)p-i (i = 0, 1, ■■■,«-!). 
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The coefficients b..(x) which appear here are continuous together with all their 
derivatives, and 

(19) b..(x) = a f .(x) (i=o, 1, •..,»— l;i=0,l, •••, m). 

The homogeneous linear differential equation of order n with solutions u.(x, p) is 

gW g[»-»] ... a 

«{•!(», p) «f— 1 ](a»,p) ••• u x {x,p) 



(20) 



For the elements of this determinant we have 



= 0. 



«ffl(«,p)-« v - E W*)^ 

4=0 



(21) 

where by (10) 

(22) \ {J0 (x)=[w.(x)yu. (x). 

Thus if we factor out of (20) 

the differential equation takes the form 

(23) % (x, P )zK + £_>, p)»W + •■• + ^(as, p)» = 0, 

where /3.(a:, p) are polynomials in 1/p. We have for the principal term of 



/U») - 



"1,11-1,0 "l.n-J.O 



^2,)i-l,0 ^"2,11-2,0 



\, tt-1, \, n-2, 



In view of (22) this last determinant may be written 

{w 1 (x)y- 1 [»!(«)]—• ••• w t (x) 1 
[w.faj)]- 1 [^ 2 (c«)] B - 2 ••• «,(*) 1 



II ««(*)■ 



[".(as)]- 1 [".(as)]-' ... w.(a,) 1 
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This is not zero at any point on (a, b) as the w t (x) are distinct and the 
u m (x) =^ by Lemma I. Therefore for \p\ = i?° 



\Pi(«>p) 



^ »>'')-^i-s^ 



But we can write (23) in the form 

(24) *M + 6 n _ t (x, p )«[-»] + ... + &„(«, /»)*=<>• 

Now the functions &..(*) and their derivatives are continuous since the to^as) 
and u t (x) are of this character. The first part of the lemma is thus proved. 
Now let j he the smallest value of j for which, for some i and x , 

From (21) and (22) we see that the principal term of 



(25) B [«.(», p)-] - A [ Ui (x, P )]=B-A [>,.(*, ,)] 
is 

(26) r3 °[% [ 5 «o( aj )-°v.( 9! )] l><( as )]']««.( a ')- 

Assume _/ = »i if possible. In each part of the difference (25) the coefficient of 

p / Wi(t)dt 

e Ja p-i» 

must then vanish, in the first part since u.(x, p) are solutions of -Z?(s) = 0; 
in the second part, by Lemma I. Therefore from (26) 

*=° 

We conclude that 

Ko( x ) = a *;,( x ) (* = 0, 1, ■■■,»! — 1), 

the «J 4 (a!) being distinct. This is a contradiction. Hence _; >• m. 

Theorem. On. a region S there exist n independent solutions 

^ + po,_ 1 (». P) ^izi + • • • + f<* (x, p)z = 
analytic in p such that if the integer m is chosen at pleasure and p is on S, 

P /""WO* rt 

z.(x,p) = u.(x,p)+eV<' -E oP - m , 
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(27) S"«(*' ') = lc^ x ' ') + e ' ^ + ' 



where 

u.(x, p)= e Ja 2-, u iA x )P 

and u ia (x) does not vanish at any point of (a, b). 

Proof. We proceed to effect a comparison of the solutions of the given dif- 
ferential equation with the solutions of the differential equation of Lemma II. 
First we write the given differential equation in the form 

( 28 > ^~n + pK-ri^, P)^A+ ■ ■ ■ + P n \(^, P)^ = p n B^A(z). 



The development of the coefficient of z^ in B — A(z) begins with a term in 
p~ m ~\ in view of (18) and (19), so that 

(29) \b k (x,p)-a h (x,p)\^B.\p\-^ 

for \p\mB . 

The general solution of a non-homogeneous linear differential equation 

< 30 ) £ + *«-:(*) |3 +---+*o(*)2/ = 4>(*0 

may be written* in terms of a set of linearly independent solutions 
y x (x), y 2 (x), • ■ ■ , y n (x) of the reduced equation in the form 

(31) y = t,c i y i (x)+ rrZy J (sH(£)'U(f)#. 

where the functions z t (x), z 2 (x), ■ ■ ■ , z n (x) are determined from the conditions 



£(a?*<">) ■*<*>- II 



(1 = 0,1, •■•,n — 2), 

(/ = »-!), 



while Cj, c 2 , • • •, c n are arbitrary constants. Any y and <f> satisfying (30) also 
satisfy (31) for some choice of c { , c 2 , • • • , c n ; and if y and <f> satisfy (31) for 
some Cj, c 2 , •••, c n , they also satisfy (30); hence (30) and (31) are entirely 
equivalent. 

*Cf., for instance, SCHLKSINGEK : Handbuch der Theorie der linearen Differentialgleichungen, 
toI. 1, p. 78. 



(J = 0,l,--.,n — 2), 

p~ n + 1 (Z = n — 1). 
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Thus treating (28) as a non-homogeneous equation we find an equivalent 
equation 

(32) z(x,p) = ±c i u i (x, P ) + r\± Ui {x, p)^!-, P )\p»B=A\z{i;,p)-\ai; 

i=l *J* L*=l J 

where the v?& are obtained from the equations 

T,u\ n (x,p)v.(x,p) = \ 

If in place of these functions v t we introduce functions u i defined by the equation 

v i (v,p)=P~ n+1 u i (x,p), 
equation (32) becomes 

n /»x f~ n ~l 

(33) z(x,p) = Y,c i u i (x,p) + pi \j^u i (x,p)u i {£,p)\B-A[z(Z,p)]d%, 
where the u. are determined from the equations 

(34) t «w (x, P )u i{ x, P )=[i wzi'X': n ~ 2h 

In order to prove the existence of a solution z k (x, p) of the character stated 
in the theorem for some definite k we make the final transformation of the 
constants c 4 in (33), 

c.= c[ (»' = l,8, ••-,*), 

( ) c t =-p J « 4 (f, f>)5lTI[a(f, p)]rff + c; (i=*+l,* + 2, •••,»)• 

This transformation is reversible. The given differential equation thus appears 
finally in the form of an equivalent integral equation 



(36) 



J ( x > P) =Z <«<(*» P) +P\ jl^x, />)«$» P) -B - 4[ss(f , p)]rff 
<=i i/« L»=i J 

+ /> ff i: ^ p)^(i, P )\b'^a[z(i, P )]rff , 



that is, the solutions of (36) and J. (s) = are the same. In the form (33) we 
could infer one and but one solution z(x, p) for a given set of c,, c 2 , • • •, c n ; in 
fact, that solution of A(z) = which satisfies the conditions 

dp n d* 

^ 7 s(«,p) = X)c.^.M i (a, p) (i=o,l f .-- t n-l), 

but a similar inference is not possible for (36). 
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For the consideration of (36) it is fundamental to establish that the functions 
<f> kl and yjr M are bounded in the following manner for | p | = IP : 

(38) (1=0,1, ••-, »— l). 

To prove this we note first that from (21) 



(39) 



l'i(x, P ) = 



n a (nB t p), 



in which v a (x, p) denotes a polynomial in 1/p. If now we substitute in (34) 

(40) u i (x, P ) = e^ WiW \(x,p), 

we find for the determination of rj^x, p) the n linear equations 

/ 41 n x-. . / ^/.... -n f" (Z = 0,1, •••,n-2), 



Tv„(x,p) Vi (x,p) = L {l=n '_i h 



The determinant of these equations is a polynomial in 1/p which we called 
/8 n (x, p), and has the principal term fi n0 (x) [see (23)] which is not zero. We 
conclude that for \p\ = JR ° 

(42) \ Vil (x,p)\S= V , \U^P)\ S V- 

Now from (39) and (40) we have 



p ft Wi(t)dt 



(43) \u¥>(x,p)u i (Z,p)\^\e 
whence 

(44) \u\H»,p)*,tf,p)\£\<'- / ?*« ) *r?\, 
if we make the restriction 

| = X (i = l, 2, ■•-, k), % = X (i — k+1, h + 2, •••, n). 

To see this one recalls the inequalities (7). Hence 



(45) 






p/J x w*(0* , 



S(n — &) 



p y '/ : w t (.t)dt 
e t) 



(fS»), 



which are in effect inequalities (38). 
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We now consider that solution z k (x, p) of (36) for which 

and we will show 

(a) that for | p | = R' one and but one such solution exists, and this is 
analytic in p ; 

(6) that the solutions z l , z 2 , ■ • •, z n thus denned fulfill the relations (27) of 
our fundamental theorem. 

The linear independence of z i , z 2 , ■ ■ • , z n is then an immediate consequence of 
their form (27), since u l , u 2 , ■ ■ •, u n are linearly independent. The demonstra- 
tion is thus completed. 

Proof of (a) . 

We know there exists one and but one solution of (33) for all sets of values 
of c { . To each set corresponds a definite transformed set c' { which we will show 
has the form 

n 

c 'i =Z7«(p)ij. 

where the 7,^(p) are analytic in p. To prove this statement we define 
Z h {x, p), (& = 1, 2, • • -, n), to be that solution of A(z) = which satisfies 
the equations 

d J d 3 

-^Z k (a,p)=j^u k (a,p) (j = o,l,...,»-l). 

If z is the solution of (32) for the set c 15 c 2 , • • •, c n , we have then 

n 

z = Hc j Z j {x,p)- 

3=1 

The Z k (x, p) are analytic in p since the coefficients in A(z) = and the 
u t (x, p) are. If this value of » is substituted in (35), we obtain the transfor- 
mation in the stated form. 

From this we see that either a unique solution for the set c\ of (46) exists, in 
which case this solution will be analytic in p, or there is a solution of (36) for 
c^ = ( i = 1 , 2, •••, n). If then we prove that the latter alternative is 
impossible for \p\ > R', we shall have proved statement (a). 

Let us now write down (36) and the equations obtained from it by differ- 
entiation, using the definitions (37). We then obtain 

*['](«., p) = ± c \ ■ «[%, P) + P Ce 9f t w * <t)dt 4> k lx, f ; p)B=A [>(f, p)]rff 
(47) "' J " 

+ P£e^ mt f kl (x, f ; p)ITi [.(f , p)] df 

(« = 0, 1, •■-, n — 1). 
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These are a set of linear integral equations in z(x, p), z^(x, p), ■ ■ •, z^"~^(x, p) 
of the Fredholm type. 

Assume if possible a solution z(x, p) to exist for c' { =0 (i = 1, 2, •••, w). 
If we write 

zW(x, p) = e'^'^O, p) (1 = 0, l, .... n-1), 

the equations (47) become 

(48) *, ( x , P ) = P f <f>jx, i-,p)g [<?, p)] rff + p f X + k! & f ; P)<7 !> (f . P)] ^ 

(J = 0, 1, ••-, n — 1) 
where 

(49) 0|>(*» p)] =Z [6,(*. P) -Oy(». />)] *,(»./>)• 

If TF=(= be the maximum of 

!«,(», /») | (?=o,i,--,ji-i) 

on (a, 6), we conclude from (29) that for \p\ = R" 

(50) \g{z{Z,p)]\^nD-W-\p\- m -K 

But in one of the equations (48), e. g., I = l x , |»,(x, p) | has the value TF at 
a: = cc, , so that 

(51) tt= L P^fa, f ; pMKI, p)]#+p f V^f; pM<& />)]<*£ 

I Ja Jb 

Applying to (51) inequalities (50) and (38), as is possible since f = x l in 
0*i, (*i» I; P) o f (38) f =a3, and in V«,( x i' £? p), we find 

(52) TT=M(6-a)r-Z> • Tripl-" 1 

which is not possible for 

\p\ ^Yn{b^a)f^n. 

Thus the set of values c\ = is seen to be impossible for | p | = R. Hence a 
unique solution for (46) exists. 

Proof of (h). 
By (46) we have 

^{x,p) = uf{x,p) + p£^uW(x,p)u i {Z,p)^B^~Aiz k {1;,p)]d1t 
(53) 

c* x r n ~\ 

+ p\ Ef^/.is^.fjtawf.p)]^ 

Writing then 
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we obtain 

(54) Ja r 

+ pj ^(^t; p)9lh(Z>p)] d %> 

where y^x, p) and g(z) are defined as in (39) and (49). This equation (54) 
is analogous to (48). 

Let W denote the maximum of 

\z M (x, p)\ (1=0,1, •■-,»»— 1) 

on (a, b). Then we find in analogy with (50) that for \p\ = H 9 

(55) \g{z h (x,p)]\^nDW\p\- m - 1 . 

If this maximum be attained for I — l^ x = x lt we find in analogy with (52) 
W^ | %! >i> P)\ + n ( b ~ a)TD- W\p\' m , 

whence, if m > and | p | is large, 

(56) " W^Q, 

since 

\Vk h (Xi,p)\ =V- 
But from (54) 

I c x 

\«ki(x, p) — y kt (x, p)\ = \p I 4>u(x, |; p)g[z k (£, p)]d£ 

] *Ja 

+ pf b *«(*.£; p)9[h(Z>p)]d£ • 

Therefore, using (38) and the inequality 

\9l*>{*>P)l\ = "D-Q\p\-~- l > 
which is a consequence of (56), we see that for a large enough | p | 

K(«. P) - %,(x, P )\ ^ n(b - a)TD- Q\p\-». 

Recalling that E is a generic notation for functions of p and other variables 
bounded for large \p\, we conclude from this at once that 



«, 



(x,p) = u k (x,p) + e>S^< l >B a -p-™, 



d , d 



/ - d , x p f x wdt)dt _-, , , 

dx Z ^ X ' P) = dx U " {X ' p) + e Bl ' P ' 






